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Four-neutrino Oscillations at SNO 

C. Pena-Garay a * 

a Inst. de Ffsica Corpuscular, C.S.I.C. - Univ. de Valencia, Spain 

We discuss the potential of the Sudbury Neutrino Observatory (SNO) to constraint the four-neutrino mixing 
schemes favoured by the results of all neutrino oscillations experiments. Our results show that some information on 
the value of cos (tf 2 3)cOS 2 (tf 2 4) can be obtained by the first SNO measurement of the CC ratio, while considerable 
improvement on the knowledge of this mixing will be achievable after the measurement of the NC/CC ratio. 



1. Introduction 

The Sudbury Neutrino Observatory is a sec- 
ond generation water Cerenkov detector using 
1000 tonnes of heavy water, D2O, as detection 
medium. SNO was designed to address the prob- 
lem of the deficit of solar neutrinos by having sen- 
sitivity to all flavours of neutrinos and not just 
to v e , allowing for a model independent test of 
the oscillation explanation of the observed deficit. 
Such sensitivity can be achievable because ener- 
getic neutrinos can interact in the D2O of SNO 
via three different reactions. Electron neutrinos 
may interact via the Charged Current (CC) reac- 
tion. All non-sterile neutrinos may also interact 
via Neutral Current (NC) and via Elastic Scat- 
tering (ES) (smaller cross section). 

The main objective of SNO is to measure the 
ratio of NC/CC events. In its first year of opera- 
tion SNO is concentrating on the measurement of 
the CC reaction rate while in a following phase, 
after the addition of MgCl2 salt to enhance the 
NC signal, it will also perform a precise measure- 
ment of the NC rate. It is clear that a cross- 
section-normalized and acceptance-corrected ra- 
tio higher than 1 would strongly indicate the os- 
cillation of v e into and/or v r . On the other 
hand a deficit on both CC and NC leading to 
a normalized NC/CC ratio 1, can only be made 
compatible with the oscillation hypothesis if v e 
oscillates in to a sterile neutrino. 
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Most of the studies of the potential of SNO 
have been performed in the framework of oscilla- 
tions between two neutrino states where v e oscil- 
lates into either an active, v e ^ or a sterile, 
v e — y v s , neutrino channel. On the other hand, 
once the possibility of a sterile neutrino is consid- 
ered, these two scenarios are only limiting cases 
of the most general mixing structure which 
permits simultaneous v e — > v s and v e — > v a oscil- 
lations. We consider those four-neutrino schemes 
favoured by considering together with the solar 
neutrino data, the results of the two additional 
evidences pointing out towards the existence of 
neutrino masses and mixing: the atmospheric 
neutrino data and the LSND results. We con- 
centrate on two SNO measurements: the first ex- 
pected result on the CC ratio and the expected 
to be most sensitive, the ratio of NC/CC. The 
measurement of other observables, such as the 
recoil energy spectrum of the CC events and the 
zenith angular dependence can provide important 
information to distinguish between the different 
allowed regions for i/ e -active oscillations but they 
are not expected to be very sensitive as discrimi- 
natory between the active and sterile oscillations. 

2. Two Neutrino Mixing: Predictions for 
SNO 

The results presented in this section have been 
obtained for the oscillation parameters in the 
presently allowed regions to the solar neutrino 
problem in the two-neutrino schemes. The re- 
sults of the global fit to the data ( Ref. ||] ) in- 
clude (i) the SK data after 1117 days of operation 
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(total number of events and day and night en- 
ergy spectra ), (ii) Gallex, GNO and SAGE data, 
(iii) the Homestake data. We used the 90% and 
99% CL regions found from the global minimum. 
The fit includes the latest standard solar model 
fluxes, BP00 model [Q. For details on the statis- 
tical analysis applied to the different observable 
we refer to Ref. ||. 

The total number of events in the CC reac- 
tion at SNO is calculated with the vd CC cross 
section computed from the corresponding differ- 
ential cross sections || folded with the finite en- 
ergy resolution function of the detector and inte- 
grated over the electron recoil energy. For defi- 
niteness,we adopted the most optimistic total en- 
ergy threshold Eth = 5 MeV. The total number 
of events in the NC reaction at SNO is obtained 
using the vd NC cross section from ||. 

In order to cancel out all energy independent 
efficiencies and normalizations we will use the ra- 
e [CC] where N^ M is the pre- 
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dieted number of events in the case of no oscilla- 
tions. The equivalent expression for the NC ratio 

-R^vc — ^ss c m = [NC] Out of those ratios one can 



compute the double ratio ^ = [NC]/[CC] for 

which the largest sources of uncertainties cancel 
out Q. As it was shown in Ref. Q, the ratio 
between the NC and CC reaction cross sections 
is extremely stable against any variations of the 
inputs of the calculations. The expected total un- 
certainties for the [CC] ratio and the [NC]/[CC] 
ratio are 6.7 % and 3.6 % respectively assuming 
5000 CC events and 1219 NC events §. 

In Fig. we show the predicted [CC] and 
[NC]/[CC] ratios for the allowed regions in the 
two flavour analysis. The dots correspond to the 
local best fit points and the error bars show the 
range of predictions for the points inside the 90 
and 99 %CL allowed regions. The mapping of the 
regions onto these bars can be easily understood 
from the behaviour of the probability for the dif- 
ferent solutions: 

(a) For oscillations into active neutrinos the 
[NC]/[CC] ratio is simply the inverse of the [CC] 
prediction. Therefore: in the SMA region smaller 
mixing angles are mapped onto higher (lower) val- 
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Figure 1. [CC] and [NC]/[CC] predictions at SNO 
for the allowed regions in the two-neutrino mix- 
ing scenarios obtained from the global analysis of 
solar neutrino data at 90 % and 99 %CL. 



ues of [CC] ([NC]/[CC]) ratio. One may notice 
that the prediction for the [CC] rate for the global 
best fit point (0.60) is larger than the measured 
rate at Super-Kamiokande. This is due to the 
nearly flat spectrum at Super-Kamiokande which 
implies that the best fit point in the global anal- 
ysis corresponds to a smaller mixing angle than 
the best fit point for the analysis of rates only; in 
the LMA region, the lower Am 2 and 9 values are 
mapped onto higher (lower) [NC]/[CC] ([CC]) ra- 
tios and viceversa; in the LOW region the higher 
(lower) [NC]/[CC] ([CC]) ratio occurs for smaller 
9 and higher Am 2 . 

(b) For the sterile case, the best fit point in 
SMA occurs at lower Am 2 than in the active case 
and this produces a higher prediction for the [CC] 
ratio (0.76). The [NC]/[CC] ratio takes an almost 
constant value very close to one (0.98 in the best 
fit point), since both numerator and denominator 
are proportional to (P v ^ Ve ). It is smaller than 
one because for the SMA solution the probability 
increases with energy in the range of detection 
at SNO and the threshold for the NC reaction is 
below the one for the CC one. 

What we see from these results is that while 
the data on [CC] can give a hint towards large or 
small mixing solutions, it will be hard to distin- 
guish active from sterile oscillations on the only 
bases of this measurement. This is not the case 
for the [NC]/[CC] ratio where both scenarios ap- 
pear nicely separated. It is not hard to fore- 
see from these results that from the [NC]/[CC] 
measurement SNO will be able to constraint the 
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additional mixings in the four neutrino scenario 
which describe the admixture of active and sterile 
oscillations. This is the main point in this paper. 

3. Allowed Four neutrino Mixing Parame- 
ters 

Together with the results from the solar neu- 
trino experiments we have two more evidences 
pointing out towards the existence of neutrino 
masses and mixing: the atmospheric neutrino 
data and the LSND results. All these experi- 
mental results can be accommodated in a single 
neutrino oscillation framework only if there are 
at least three different scales of neutrino mass- 
squared differences. The simplest case of three in- 
dependent mass-squared differences requires the 
existence of a light sterile neutrino, i.e. one whose 
interaction with standard model particles is much 
weaker than the SM weak interaction, so it does 
not affect the invisible Z decay width, precisely 
measured at LEP. 

There are six possible four-neutrino schemes 
that can accomodate all these evidences. They 
can be divided in two classes: 3+1 and 2+2. In 
the 3+1 schemes there is a group of three neutrino 
masses separated from an isolated mass by a gap 
of the order of leV which gives the mass-squared 
difference responsible for the short-baseline oscil- 
lations observed in the LSND experiment. In 2+2 
schemes there are two pairs of close masses sep- 
arated by the LSND gap. 3+1 schemes are dis- 
favoured by experimental data with respect to the 
2+2 schemes but they are still marginally allowed. 

As discussed in Ref. ||, for any of these four- 
neutrino schemes, either 2+2 or 3+1, only four 
mixing angles are relevant in the study of solar 
neutrino oscillations The survival of so- 

lar z/ e 's mainly depends on the mixing angle #12, 
whereas the mixing angles $23 and $24 determine 
the relative amount of transitions into sterile v„ 
or active v a , this last one being a combination of 

and v T controlled by the mixing angle #34. 
and v T cannot be distinguished in solar neutrino 
experiments, because their matter potential and 
their interaction in the detectors are equal, due 
only to NC weak interactions. As a consequence 
the active/sterile ratio and the survival probabil- 



ity for solar neutrino oscillations do not depend 
on the mixing angle $34, and depend on the mix- 
ing angles $23 #24 only through the combination 
cos $23 cos #24 ■ For further details see Ref. @,| . 

In the general case of simultaneous v e — > v s 
and v e — > v a oscillations the corresponding prob- 
abilities are given by |l],|| 

Pv e —>v s = c 23 c 24 (X — Pv a — >i/ e ) j (1) 

P^u. = (1 - C2VI4) (1 - P v .^ v .) . (2) 

where Py K -n> e takes the standard two-neutrino 
oscillation form for Am^ and 612 but com- 
puted with the modified matter potential A = 
Acc + c^s^Anc- Thus the analysis of 
the solar neutrino data in the four-neutrino 
mixing schemes is equivalent to the two- 
neutrino analysis but taking into account that 
the parameter space is now three-dimensional 
(Amf 2 , tan 2 $12, cos 2 #23 cos 2 $24)- We want to 
stress that, although originally this derivation 
was performed in the framework of the 2+2 
schemes pL§], it is equally valid for the 3+1 ones 

i- 

The allowed regions in the three-parameter 
space for the global combination of observables 
are shown in Ref. 0. The global minimum 
used in the construction of the regions lays in 
the LMA region and for pure z/ e -active oscilla- 
tions, c| 3 c| 4 = 0. The SMA region is always a 
valid solution for any value of c\ z c\^ at 95% CL 
( 0.1K c% z c\ A <0.31 allowed at 90% CL). On 
the other hand, the LMA and LOW-QVO solu- 
tions disappear for increasing values of the mixing 
c| 3 ci 4 (0.72 and 0.76 at 99% CL , respectively). 

4. Expected Rates at SNO in Four- 
Neutrino Schemes 

In this section, we present the predictions for 
the CC ratio and for the NC/CC ratio in the four- 
neutrino scenario previously described. This sce- 
nario contains as limiting cases the pure ^.-active 
and i/ e -sterile neutrino oscillations. 

In Figs. ||— [| we show the results for the pre- 
dicted [CC] ratio and [NC]/[CC] ratio for the dif- 
ferent allowed regions (SMA, LMA, LOW-QVO) 
at 90 and 99 %CL as a function of c\ z c\±. The 
general behaviour of the dependence of the pre- 
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dieted ratios with c^c^ can be easily understood 
using the following simplified expressions: 

[CC] ~ P u ^ Ve , (3) 
[NC] \-c\A(\-P v ^) 

[cc] ~ p Ve ^ Ve • 1 ' 

From Eq. (||), the only dependence of [CC] on 
c 23 c 24 1S due to the modification of the mat- 
ter potential entering in the evolution equation 
and it is very weak. The dependence of the al- 
lowed range of the [CC] ratio with c^c^ dis- 
played in the figures arises mainly from the vari- 
ation of the size of the allowed regions. Alter- 
natively following Eq. (Q) we find a stronger lin- 
ear dependence of [NC]/[CC] on c^c^ with slope 
(-1 + P Ue ^u e )/Pu e ^u e ~ 1 - 1/[CC] and inter- 
cept l/P„ e ^ Ve ~ 1/[CC]. This simple description 
is able to reproduce the main features of our nu- 
merical calculations as can be seen in the figures. 

Figs. |2|-^ contain the main quantitative result 
of our analysis in the four- neutrino mixing sce- 
nario. From each of them it is possible to in- 
fer the allowed range of the active-sterile admix- 
ture, C2 3 c| 4 , compatible, within the expected un- 
certainty, with a given SNO measurement of the 
ratios. Also, comparing the allowed ranges for the 
different solutions one can study the potential of 
these measurements as discriminatory among the 
three presently allowed regions. Of course, both 
issues are not independent as we have no "a pri- 
ory" knowledge of which is the right solution and 
both must be discussed simultaneously. In order 
to do so we pass to describe and compare in detail 
the predictions in the different regions. 

The results for the SMA solution are shown 
in Fig. |a and |.b for [CC] and [NC]/[CC] ra- 
tios respectively. First we notice that we find 
a small region allowed at 90 %CL only for a 
non-vanishing admixture of active and sterile os- 
cillations as mentioned before. In this region 
[CC]~0.65-0.73 and [NC]/[CC]~ 1.3-1.4. The 
predictions at 99% range from [CC]~ 0.4-0.9 
([NC]/[CC]~l.l-2.5) for pure i/ e -active scenario 
to to [CC]- 0.59-0.85 (fNC]/[CC]~0. 96-0.98) for 
pure v e -sterile oscillations. Thus if SNO ob- 
serves a ratio [CC]<0.58 the value of c^c^ can 
be constrained to be smaller than 1 disfavouring 
pure ^ e -sterile oscillations. On the contrary a 




COS 2 1?23COS 2 "# 2 4 COS 2 T?23COS 2 "l9 2 4 



Figure 2. [CC] and [NC]/[CC] predictions at 
SNO as for the SMA region in the four-neutrino 
scenario obtained from the global analysis of so- 
lar neutrino data at 90 % (lighter) and 99 %CL 
(darker). The dotted line corresponds to the pre- 
diction in the case of no oscillations. 



measurement of [CC]>0.68 will immediately hint 
towards the SMA solution but will not provide 
any information on the active-sterile admixture. 
Also, one must notice, that such value, although 
allowed by the present global statistical analy- 
sis at 99 %CL, will imply a strong disagreement 
with the total rate event rate observed at Super- 
Kamiokande. 

As seen in Fig. |.b the [NC]/[CC] ratio is 
more sensitive to the active-sterile admixture. To 
guide the eye, in the figures for the [NC]/[CC] 
ratio we plot a dotted line for the prediction in 
the case of no oscillation [NC]/[CC]=1. For any 
of the solutions, the allowed range for this ra- 
tio shows as general behaviour a decreasing with 
c 23 c 24 due to two effects: (i) the allowed regions 
become smaller and (ii) the prediction decreases 
when more sterile neutrino is involved in the oscil- 
lations as described in Eq. (||). The measurement 
of higher values of this ratio will favour the four- 
neutrino scenario with larger component of in- 
active oscillations. On the other hand a measure- 
ment of [NC]/[CC]~ 1, will push the oscillation 
hypothesis towards the pure v e -sterile oscillation 
scenario. This case will be harder to differentiate 
from the non-oscillation scenario. We find that 
with the expected sensitivity the parameter c| 3 c| 4 
is constrained to be above 0.44 at 99% CL and 
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that the pure v e -&ctive oscillations in the SMA 
region are compatible with [NC]/[CC] = 1 only 
at ~ 5(7. The predictions for oscillation param- 




Figure 3. Same as Fig. | for the LMA region. 



etcrs in the LMA region are shown in Figs. |[a 
and |.b for [CC] and [NC]/[CC] ratios respec- 
tively. The predictions at 99% vary in the range 
[CC]- 0.18-0.62 and [NC]/[CC]~ 1.4-5.6. The 
first thing we notice by comparing Fig. |^.a with 
Fig. ||.a and Fig. ||.a is that the most discrimina- 
tory scenario for the [CC] rate results if SNO finds 
a small value [CC]— 0.25. This would signifi- 
cantly hint towards the LMA solution to the solar 
neutrino problem and towards and i/ e -active os- 
cillation scenario. First, it is well separated from 
the predictions for the SMA and LOW regions. 
Second, it will include as a bonus a small but mea- 
surable day-night asymmetry. And third it will 
constrain the c| 3 c| 4 to a small value (~ 0.2). On 
the contrary the less discriminatory scenario will 
be a measurement 0.4<[CC]<0.6 where the pre- 
diction would be compatible with both SMA and 
LOW-QVO solutions and no improvement on our 
knowledge of the four-neutrino schemes is possi- 
ble. The [NC]/[CC] ratio can definitively improve 
the discrimination between the different scenarios 
provided its measurement lays in the upper range. 
For instance a measurement of [NC]/[CC]~ 4 (± 
0.7 at 5cr) will be conclusive for selecting LMA as 
the solution to the SNP and will imply and upper 
bound on c\^c\ A < 0.3. The predictions for the 
LOW-QVO region lay between the ones for SMA 



hiO.E 
O.f 
0.1 

o.t 

0.1 
0.( 

o.: 
o.i 

0.1 



0.2 0.4 0.6 0.8 




0.2 0.4 0.6 0.8 • 



Figure 4. Same as Fig. | for the LOW region. 



and LMA as displayed in Fig. [| and therefore 
they are more difficult to discriminate. The pre- 
dictions at 99% vary in the range [CC]— 0.3-0.68 
and [NC]/[CC]~ 1.2-3.4. As a consequence we 
see that a low [CC] ratio but still within the 99% 
CL range allowed for this region, 0.3<[CC]<0.4, 
will constrain significantly the c| 3 c| 4 parameter 
compatible with this solution but it will not be 
distinguishable from the LMA solution unless the 
measured [CC]<0.3. As mentioned above, the 
[NC]/[CC] ratio will be able to differentiate the 
LMA and LOW-QVO solutions if not in the range 
[1.5,3]. One should also notice that for the upper 
part of this range a positive measurement of the 
day-night asymmetry and the zenith dependence 
will point towards the higher Am 2 of the LOW 
region as the solution. 
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